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ABSTRACT 

The recent detection of a transient absorption feature in the X-ray prompt emission 
of GRB 990705 showed the importance of such observations in the understanding of 
gamma-ray bursts and their progenitors. We investigate the time dependence of pho- 
toionization edges during the prompt emission of bursts in different environments. We 
show that their variabihty can be used to infer the density and geometry of the sur- 
rounding medium, giving important clues to unveil the nature of the burst progenitor. 

Key Virords: Gamma-rays: bursts — X-rays: general — X-rays: ISM — line: forma- 
tion 



1 INTRODUCTION 

The nature of the progenitors of gamma-ray bursts (GRBs) 
is still mysterious, since the fireball producing the gamma- 
ray emission does not carry any information about the pro- 
genitor generating it (see, e.g., Piran 1999). Recent ob- 
servations (e.g., KuLkarni et al. 1999) suggest the associa- 
tion of GRBs with the final stages of evolution of massive 
stars (Woosley 1993, Paczyhski 1998; MacFadyen & Woosley 
1999), but many crucial points still remain controversial. 
The possible evidence of rebrightening shown by the optical 
afterglows of several bursts ~ 30 days after the burst explo- 
sion (Bloom et al. 1999; Reichart 1999; Galama et al. 2000) 
suggests a simultaneous explosion of the burst with a Type 
Ic supernova. On the other hand, the detection of iron ab- 
sorption and emission features in the X-ray emission of the 
bursts and afterglows (Amati et al. 2000; Piro et al. 2000) 
can be explained with a two-step explosion, in which the 
burst onset follows a supernova explosion by several months 
(Lazzati et al. 1999; Vietri et al. 2001; Lazzati et al. 2001). 

These different scenarios can be distinguished through 
the analysis of the radial distribution of the interstellar 
medium surrounding the explosion site (Ghisellini et al. 
1999). If the burst explodes simultaneously with the su- 
pernova, the photons propagate through the pre-explosion 
stellar wind, while in the two-step scenario a high density 
metal enriched supernova remnant is expected to surround 
the burst explosion site (Lazzati et al. 1999). 

In this paper we study the variability properties of pho- 
toionization edges in the first several tens of seconds from 
the burst onset, for the different radial distributions of the 



density profile characterizing these two scenarios. Temporal 
variations of the opacities are indeed expected as a result 
of the gradual photoionization of the medium by the X-ray 
prompt emission of the burst. A similar variability effect on 
absorption lines had been discussed by Perna & Loeb (1998) 
and, in the X-rays, by Bottcher et al. (1999). In these pa- 
pers, however, dependences on the radial profile of the den- 
sity had not been considered. 

The paper is organized as follows: in §2, we present 
some analytic approximations for the time-dependent opac- 
ities which hold for the optically thin case, discussed in §3. In 
§4, we present the results of the general (numerical) solution 
to the problem, and compare it with the analytic approxi- 
mations derived in §2. In §5, we discuss possible effects of 
the fireball expansion in the absorbing medium. Finally, our 
results are summarized in §6. 



2 ANALYTIC THEORY 

Consider a distribution n(r) of absorbers, characterized by a 
photoionization cross-section g. In the case of an impulsive 
illumination by a strong ionizing fiux the density n(r) be- 
comes time dependent, and so does the absorption opacity: 

/>oo 

r(t) = o- / n{r,t)dr (1) 
Jo 

To solve Eq. 0, we must determine the time-dependent den- 
sity n{r, t). We use a particular time reference: at each radius 
r from the photon source, the time is set to be zero when the 
first ionizing photon passes by. By saying that the ionization 
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time at radius f is t, we then mean that t is the time interval 
between the first ionizing photon crossing the radius f and 
the complete ionization of all absorbing ions located at the 
same radius. The advantage of this time reference is that the 
travel time of photons is null, and the time t is also the time 
measured by an observer at infinity. The same conceptual 
scheme applies in §5 for the radius of the fireball. 

Consider now a medium in which recombination is not 
efficient and let us neglect the change of the resonant fre- 
quency as an atom is progressively stripped off its electrons 
from the ground state to complete ionization. The first con- 
dition is easily fulfilled in GRBs, given their high ionization 
flux during the prompt phase (exceptions may however be 
possible, see Lazzati et al. 2001). In this case the density 
n(r, t) remains constant until all electrons are stripped. The 
time of complete stripping depends on the distance r of an 
atom from the photon source. We assume that, at a given 
radius f , n(f , t) changes abruptly in time: 



n(f,t) = n(f)X[o,zti„(r)]W 



(2) 



where X[a,b][^) ~ 1 for a < x < b and elsewhere, tion is the 
ionization time of a single electron and Z is the "efficient 
number" of electrons of the relative element. Z is in general 
smaller than the atomic number Z of the element, since 
other processes can contribute to ionization. These include 
coUisional ionization, the Auger effect and the ionization of 
the external electrons by the UV flux. The actual value of 
Z can be obtained through numerical simulations and is not 
universal, being dependent on the ionizing UV continuum 
and the density of the absorbing medium (see § 4). The 
ionization time of a single electron is given by: 



^ion 



h V 



(3) 



where Fiv) is the flux of ionizing radiation at a given radius, 
taking into account the effect of absorption from material at 
lower radii. Eq. ^ modifies as: 



tion(r-) 



47rr 



47rr^ 

foo Liv) -ct(i') i ' n(p.t)dp / ^ J 

/ -r^ e Jo a[v)dv 



(4) 



where L(y) is the luminosity of the ionizing continuum and 
is the threshold photoionization frequency. Note that the 
ionization time, needed to compute the density of absorbers, 
depends on the density itself in a non-linear way. For this 
reason Eq. W does not have a general analytical solution. 



3 OPTICALLY THIN MEDIA 

If the absorbing medium is optically thin, the flux at a given 
radius does not depend (at least as a first approximation) 
on the absorption taking place at smaller radii. In this case 
Eq. U simplifies to: 



tioii(r-) 



47rr 



\^h^ylo{v)dv 



(5) 
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Figure 1. Opacity vs. time for a uniform cloud of absorbers. Ini- 
tially opaque, intermediate and thin clouds are considered, from 
top to bottom (see text for the parameter set). The dashed lines 
correspond to the analytic approximations of Eq. |ll| while the 
symbols are obtained with the numerical simulations described 
in §4. As expected, the analytic approximation holds in the op- 
tically thin regime, but fails to reproduce accurately the drop of 
the opacity in the optically thick cases (top line and diamonds). 



By approximating the shape of a photoionization edge as 
(t(!/) = (To {iy/uo)~^, and considering a power-law spectrum 
for the ionizing fiux [L{u) = Lq (u/i/q)^"], Eq. ^ becomes: 



tion{r) = (3 + a) /i 



Lo (To 



Combining Eq. ^ with Eq. |l] and Eq. |^ we obtain: 



-{V-.t) = (To 



1^0 



n(r) dr 



(6) 



(7) 



where Rii(t), the radius of the ionization front at time t, is 
obtained by inverting Eq. ^: 

1/2 



Lo (To t 



(8) 



^47r/iZ(3 + a), 

Assuming an absorbing medium whose density decreases as 
a power law with distance, [n{r) — uq {r / Ro)~^], we obtain 



T{i',t) = (To ( — ) no -Rq 



r dr 



(9) 



max[J?o )-Rif (*)] 



where Rm and Rq are the maximum and minimum radii 
of the density profile. This avoids unphysical divergences at 
very large and very small radii. With this assumptions Eq. ^ 
becomes: 



T{u,t) = (Tonoi?o f~~) ^ 



In 



f max[flM,-Rif(t)l 1 
I max[flo,flif(t)l J 



max[flM-Rif(t)l^~^^-max[flo.flif(t)]^-3 I 



(10) 



Two particularly interesting cases are a uniform distribution 
up to a given radius Rm and a rapidly decreasing (/3 > 1) 
distribution extending to infinity. 
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3.1 Uniform density 

In the case of a uniform density of absorbers {(3 — 0), Eq. |lC 
becomes particularly simple: 



aono 




Kit < Ro 

Ro < Rif < Rm 

Rii > Rm 



where A is given by: 



A 



A-K h Z {2, + a) 



(11) 



(12) 



Note that the case of uniform density can correspond not 
only to the case of a burst surrounded by a standard inter- 
stellar medium, but would also be a reasonable approxima- 
tion for the case of a young supernova remnant (see Lazzati 
et al. 1999). 

The opacity is constant at the beginning, then de- 
cays becoming rapidly negligible. Examples of this behav- 
ior are shown in Fig. |l| (dashed lines) for three different 
initial conditions: an opaque (r(0) ~ 13) cloud, an in- 
termediate (t(0) ~ 1.3) cloud and a thin (r(0) ~ 0.13) 
cloud. We have considered the photoionization edge of iron, 
with a resonance cross section ao ~ 2 x 10~'^° cm^ and a 
threshold frequency hi/ = 9.28 keV. The parameterization 
adopted for the ionizing flux and the ambient medium is: 
Lo = 1.75 X 10^^ erg s~^ Hz"^, q = (similarly to the spec- 
trum of GRB 990705, Amati et al. 2000), Ro = lO" cm, 
^iM = 1 pc and no = 200, 20, 2 cm~"^, respectively. These 
conditions are quite extreme, since they correspond to a to- 
tal mass of the absorbing iron of ~ 1000, 100 and 10 Mq, 
respectively. On the other hand, a r ~ 1 feature lasting for 
~ 10 seconds has indeed been detected (Amati et al. 2000), 
and very large iron masses are required, in some scenarios, 
to explain this feature (see Lazzati et al. 2001 and references 
therein for a more complete discussion). 

We have considered in the simulation a constant ioniz- 
ing luminosity. This is not generally true for GRBs, where 
fluctuations on timescales of fractions of seconds are usually 
observed. This would cause small scale fluctuations, with 
the same timescale, on the time evolution of the opacities, 
but the general trend is unaffected. Moreover, the burst may 
turn off before the complete ionization of the surrounding 
medium. In this case the results of this paper would hold 
only for times smaller than the turn off time of the GRB. 

By comparison with the numerical simulations pre- 
sented in §4 (where the approximations of §2 are released), 
the "efficient number" of electrons of iron in these condi- 
tions turns out to be Zpe ~ 7. The dashed lines have been 
computed adopting this value. 

In Fig. ^, we show the time of complete ionization of 
the iron ions as a function of the distance from the source of 
the ionizing photons. In the analytic approximation of Eq. ^ 
this time does not depend on the initial opacity and density 
distribution of the medium. 



3.2 Decaying density profiles 

Let us consider now the case of a density profile decreasing 
with distance from the burst source. This can correspond, 
for example, to a pre-explosion stellar wind if the bursts are 
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Figure 2. Time of complete ionization vs. radius for the same 
geometrical setup of Fig. ^. The dashed line shows the analytic 
approximation of Eq. H. The symbols refer to the same cases of 
Fig. ^ Since in the simulations the approximation of the step 
function (Eq. ^ is not used, the symbols are relative to the time 
at which the fraction of iron completely ionized is a half of the 
total. The thin solid lines in the upper left corner of the figure 
show the time at which the fireball crosses the radius ijjf (see §5). 
Since this time is always much larger than the time of ionization, 
the hydrodynamics of the fireball does not infiucnce the evolution 
of the opacity. 



associated with the death of massive stars. In this case we 
would have /3 = 2 (ChevaUer & Li, 1999). 

If the density of the absorbing medium decreases with 
distance sufficiently rapidly (/? > 1), the maximum radius 
can be set to infinity and Eq. ^ becomes: 

I (At)^ Ru>Ro 



(13) 



where A is given in Eq. |l2| After a short initial constant 
phase, the opacity decreases in time as a power-law whose 
slope depends only on the index (3 of the density profile. 
The time evolution of the opacity in a wind environment is 
shown in Fig. ^. We used Ro — 10^^ cm and initial densities 
no = lO'', 10* and lO'^cm"^ (corresponding to r(0) ~ 200, 
20 and 2). The parameters describing the ionization flux 
are the same used for Fig. |l] In Fig. we show the time 
evolution of the radius of complete ionization, for the same 
parameters used for Fig. [l. 



4 NUMERICAL RESULTS 

To solve the problem in its full generality, we remove the ap- 
proximations made in § 2. We start the simulation at t = 0, 
and let the burst photons propagate. In propagating from a 
point at position r to another point at position r + Ar, the 
ioinzing flux is reduced according to: 



F4r + Ar,t + At) = F^(r,t)exp[-AT„{r,t)]- 



'(r + Ar)2 -^-^"^^ 

The optical depth due to photoabsorption within the dis- 
tance Ar is then given by: 
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Figure 3. Opacity vs. time for a wind environment. Initially 
opaque, intermediate and thin winds are considered, from top to 
bottom (see text for the parameter set). The dashed lines corre- 
spond to the analytic approximations of Eq. while the symbols 
are obtained with the numerical simulations described in §4. As 
expected, the analytic approximation holds for t < 1. 



Ar,(r,t) = Ar^nj(r,t)a,(i/) 



(15) 



The photoionization cross sections are taken from Reilman 
& Manson (1979). The ionic concentrations are determined 
by solving the system of equations: 

drij (r, t) 



dt 



= qj-2nj-2 + qj-iUj-i + Cj-iUj-me 
— {qj + CjUe + ajne)nj + aj+iUj+iTie . (16) 



The qj and Cj are respectively the photoionization and col- 
lisional ionization coefficients of ion j, while Uj is the re- 
combination coefficient. Note that qj-2 refers to inner shell 
photoionization followed by Auger ionization. The coUisional 
ionization rates are calculated according to Younger (1981). 
We compute the terms due to photoionization by integrating 
Li/CTi/ numerically. The recombination rates are given by the 
sum of the radiative and dielectronic recombination rates. 
The code uses routines developed by Raymond (1979; see 
also Perna, Raymond & Loeb 2000 for a similar application 
and further details). After updating the ionization fractions 
at each time step At, the optical depth T^{t) — J drA.T,^{r, t) 
is recomputed. 

A comparison between the numerical and analytical re- 
sults is shown in Fig. |l| |^, ^ and ^. In all these figures, the 
analytical approximation is drawn as a dashed line, while the 
numerical results are overlaid as symbols (circles, diamonds 
and triangles). The introduction of recombination and col- 
lisional ionization does not change significantly the results 
with respect to the analytic case, where only photoioniza- 
tion is considered. As expected, corrections must be applied 
in the optically thick regime. 

In Fig. 1^, the opacity is dominated by material at large 
radii. The analytic approximation is valid when the initial 
opacity is small or intermediate (r(0) ^ 1), but underesti- 
mates the time of decline for an initially opaque cloud. For 
both the initially thin and thick cases, the analytic approx- 
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Figure 4. Same as Fig. Elbut for the wind environment of Fig. 



imation declines faster than the numerical result. This is 
due to the approximation of Eq. which severely underes- 
timates the absorption at a given radius for t > iion('')- 

Fig. 1^ shows the time of complete ionization for the 
same parameter set of Fig. ^. In the numerical simulations 
the time of complete ionization is not a well defined quantity, 
since at each radius and at any time there is always the pos- 
sibility (albeit very small) to have ions with bounded elec- 
trons. To properly compare with the time defined in Eq. ^ 
which is plotted as a dashed line, the symbols (circles, dia- 
monds and triangles) show the time at which about 50 % of 
the iron is completely stripped (FeXXVII). The figure shows 
that the analytic approximation reproduces accurately the 
numeric results, with small deviations at large radii for the 
initially opaque case. 

In the case of a wind environment (Fig. ^), the opacity is 
dominated by material at small radii. The analytic approxi- 
mation describes accurately the temporal decay of the opac- 
ity in the optically thin and intermediate cases (r(0) 1). 
For an initially optically thick medium, the numerical solu- 
tion deviates at small times, but collapses onto the analytic 
approximation when the medium becomes optically thin. 
Again, Zyc ~ 7 is adequate. The analytic time of complete 
ionization for a wind environment is compared with the re- 
sults of numerical simulations in Fig. ^ Again, the analytic 
approximation is adequately accurate at all times but for 
the initially optically thick case, in which deviations are ob- 
served. 



5 FIREBALL DYNAMICS 

For all our calculations we have assumed that the ionization 
flux comes from a point source. This is not strictly true, 
since in the fireball model the prompt radiation is produced 
by the expanding shell. To check whether our assumption is 
reasonable we must demonstrate that the fireball radius is 
smaller than the ionization radius at all times. 

The radius of the relativistic fireball expanding in 
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an inhomogeneous external medium is given by (see, e.g.. 
Meszaros, Rees & Wijers, 1998): 



R{b{t) =Ris + crl 



3-|3 



t < tES 
t > tES 



(17) 



where To is the asymptotic Lorentz factor of the fireball and 
Ris ~ 10^^ cm is the internal shock radius, i.e. the radius 
at which the first ionizing photons are produced. The time 
tES is the onset time of external shocks, where the fireball 
expansion starts to be slowed down by the interaction with 
the external medium. This is given by (Meszaros, Rees & 
Wijers 1998): 



tES 



{3-f3)E 



47r rg nip c2n{H,o} 



Rr, 



3^ 



(18) 



where E is the total fireball energy, nip the proton mass, 
and a particle distribution nH{R) = n{H.o}{r / Rx))~^ has 
been assumed. Note that the distribution of particles in the 
external medium follows the distribution of the absorbers 
given above. The normalization are instead different. For 
solar iron abundance we have no — 4.68 x 10~^n{£f,o} (An- 
ders & Grevesse 1989). 

In Fig. ^ and ^ the time at which the fireball crosses 
the radius _Rif is plotted with thin solid lines. We assumed 
a fireball with E = lO^'^ erg and Fq = 100, and an iron 
abundance ten times solar. It can be seen that the fireball 
crossing time is always much larger than the time of ioniza- 
tion. We conclude that, for a reasonable set of parameters, 
the assumption of decoupling between the ionization and the 
hydrodynamics is valid. 



6 CONCLUSIONS 

We have calculated the time dependence of the absorption 
opacities due to the material in the surroundings of gamma- 
ray bursts. We have been able to qualitatively reproduce 
the numerical results with analitic expressions using some 
approximations appropriate for the optically thin regime. 

In the optically thick regime, corrections must be ap- 
plied. In the case of a uniform medium, the analytic expres- 
sions fail to reproduce exactly the lifetime of the absorption 
feature (i.e. the timescale of constant optical depth), while 
in the wind case the analytic expressions are in good agree- 
ment with the numerical results at all but the very early 
times. 

An important issue is the observability of these features. 
On one hand, we require the density of the surrounding 
medium to be high, in order to have the highest possible 
opacity. On the other hand (Lazzati et al. 2001), the Thom- 
son scattering optical depth, tt, of the absorbing material 
must be less than one, in order to maintain the fiickering be- 
havior of the 7-ray lightcurve. For a given Thomson opacity, 
the iron opacity is given by rpe ~ LSttAfb, where Aeb is the 
iron abundance in solar units. Iron enriched media are hence 
necessary to observe a deep iron feature in a tt ~ 0.1 cloud. 
Enriched media around GRBs have indeed been observed 
(Amati et al. 2000; Lazzati et al. 2001). 

For the wind environment, the absorption feature has 
an extremely short lifetime, not detectable with the present 
and near future detectors and instruments. It then appears 



that a stellar wind, even if strong and appropriate to a very 
massive progenitor, is not able to imprint a detectable ab- 
sorption feature in the X-ray spectrum of bursts. The reason 
is the immediate ionization of the absorbing material. How- 
ever, there is still the possibility that the wind of the pro- 
genitor, especially in the last phases, is unsteady and inter- 
mittent, therefore populating the burst environment with a 
density profile different from R~^ and corresponding to more 
mass at larger radii than the case we assumed. If this case 
can be described by a larger value of Ro, then the lifetime 
of the absorption feature should scale as R^ and detectable 
features may appear if Rq ~ 10^^ cm and/or with a weaker 
ionization flux, i.e. in X-ray poor bursts. If, beyond Rq, the 
density profile retains a power law profile, then also the ab- 
sorption optical depth decreases in time as a power law, and 
this discriminates the wind case from the uniform density 
scenario, where the fall off of the absorption feature should 
be more abrupt. 

The required time resolved X-ray spectroscopic obser- 
vations are beyond the capabilities of present instruments 
and satellites, which can provide only a time integrated 
measurement of the opacity on timescales of 5 10 sec- 
onds (cfr. Amati et al. 2000). The Swift satellite, however, 
will be able to slew on target and to start the observations 
with the X-ray telescope in a few seconds (in particularly 
favorable cases), and will be sensitive enough to measure 
time-resolved opacities (Perna et al., in preparation). 
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